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Abstract 

^ ' 

*l2 ' Although energy levels are often given by solutions of the radial equa- 

tion such that u(0) is non zero, and hence by first-order singular func- 

rT-l tions which are not eigenfunctions of H, the latter is always considered 

as the only operator that gives energy levels. Vibrational levels of di- 
atomic molecules are a usual example. We show that the operator which 
has singular eigenfunctions, or pseudofunctions, that give energy levels, 
is the operator whose action on pseudofunctions amounts to the embed- 
ding in the distributions of R 3 of their Hamiltonian in R 3 /{0}. When 
Cq ■ its eigenfunctions are regular, this operator amounts to H. Energy levels, 

VQ \ which are given by eigenfunctions of H when u(0) is zero, are thus given 

in any case by eigenfunctions of this operator, which is an extension of 
the Schrodinger operator, but not of the Hamiltonian. 



1 Introduction. Energy levels and radial equa- 
tion. 

In a central potential V(r) energy levels £ n are eigenvalues of the Hamiltonian 
that belong to normalized eigenfunctions tf> n = [w n (j)/r]Y^ [9, (p) which behave 
at the origin like r e Y^ (9, <p), and hence eigenvalues of the radial equation that 
belong to normalized eigensolutions w n (r) which behave at the origin like r e+1 , 



h 2 d 2 e(£+l)H 2 



2m dr 2 2mr 2 



•V(r) 



W n (r) = £ n W n (V) (1) 



As physical potentials V(r) are not exactly known, one must substitute 
theoretical and then approximate potentials V{r) for which the Hamiltonian 
has a discrete set of normalizable eigenfunctions X P„= \u n (r)/r]Y^ (9,<p), and 
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hence the radial equation a discrete set of normalizable analytic eigensolutions 
u n (r)which behave at the origin like 



r i+i 



2m dr 2 



i)h 2 



2mr 2 



+ V{r) 



(r) = E n u n (r) 



(2) 



One assumes indeed that if V(r) s=s V(r), "J/n will be similar to ip n , u n (r) to 
w n (r), and E n ss £ n . The solutions of the radial equations, which are written 

, are given in R 3 /{0} by both roots 



fc>O u fc 



dJlr k 



in the form of series, u n (r) = r Yl 
A= £ + 1, and A= -^, but in R 3 for £ > only by the root A= £ + 1 [1]. These 
solutions are normalizable and give, when substituted in \IJ n , eigenfunctions of 
H, or solutions of the Schrodinger equation. For £ = the radial equations 



2to dr 2 



w„ (r) = £„w„ (r) 



(3) 



2?ri dr 2 



+ V(r) 



U„ (r) = B„H n (r) 



(4) 



have solutions in R 3 given by both roots, that is, A = 1, and A= 0. Both roots 
give normalizable solutions, but only the solutions given by the root A = 1 give, 
when substituted in \l/„, solutions of the Schrodinger equation [1]. Insofar as 
energy levels £ n are eigenvalues of the Hamiltonian, the solutions w n (r) of (3) 
behave at the origin like r. In order for ^„ to be eigenfunctions of H, one must 
substitute potentials V(r) for which the following conditions must be satisfied: 
i) V(r) must be close to V (r), at least in the neighborhood of the minimum; 
ii) Eq.(4) must have a discrete set u n (r) of normalizable analytic solutions ; 
in) these solutions must behave at the origin like r. 

Now, explicit analytic solutions of (4) can be derived only for a limited 
number of potential energy functions V(r). When we demand that the first 
two conditions be satisfied, this number is much more limited. If we moreover 
demand that these solutions behave at the origin like r, the problem cannot 
necessarily be solved. This is why one substitutes potentials V(r) close to V (r), 
for which (4) has a discrete set of normalizable solutions regardless of their value 
at the origin. One intuitively expects indeed that if V (r) ss V (r) and if (4) has 
a discrete set of normalizable solutions u„ (r), moreover demanding that they 
behave at the origin like r is a supplementary condition that is not required to 
obtain eigenvalues E n of (4) close to the eigenvalues £ n of (3). This is shown 
by the following example. 

2 Vibrational levels of diatomic molecules. 

The vibrational levels of diatomic molecules are of the form 



£, 



■ V m + (n + \)hw + C 2 (n + i) 2 + C 3 {n + i) 3 + 



0,1,2... 



In the Born-Oppenheimer approximation the potential energy of interaction 
between the two nuclei of a diatomic molecule is a central potential V(r), at- 
tractive at large distances, repulsive at short distances, with a minimum at r = 
r m . By expanding the potential in powers of r - r m in the neighborhood of the 
minimum by neglecting the terms of order > 3, it gives 

V(r)^hmcu 2 (r-r m ) 2 -V m u = [V" {r m ) /m] 1 ' 2 (5) 

This parabolic potential is used as a first approximation, which holds in 
the neighborhood of the minimum, and then gives the first levels. The radial 
equation for £ = involving this potential, 



1 2/ \2 T ^ 

"7T— ~r^ + i moJ { r - r m) - Vrr, 

2m dr z 



u n (r) = E n u n (r) (6) 

has the square integrable solutions and the corresponding eigenvalues 



«n (r) = N n exp [-J^ 2 (r - r m ) 2 ] H n [f3(r - r m )} /3 2 = muj/h 

E n = (n + %)hu-V m n = 0,l,2... (7) 

where H n are Hermite polynomials of order n and N n normalization con- 
stants. Experimental data show that these eigenvalues give estimates of the 
first vibrational levels of most molecules. Now, except for the case where f3r m 
coincides with a zero of H n , 

u n (0) = N n exp(-hl3 2 r m 2 )H n (/3r m ) £ (8) 

Analytical forms for V(r) closer to the realistic physical potential have been 
proposed, the most frequently used is the Morse potential [2] 

V M (r) = V m {exp [-2a(r - r m )} - 2exp [-a(r - r m )}} (9) 

For £ = the square-integrable solutions and the eigenvalues of the radial 
equation are [2] 

u n (r) = exp[-z(r)/2] [z{r)f 2 L b n+b [z(r)} 
z = 2dexp [— a(r — r m )} d = (2mV m ) :2 (ah)~ b = 2d — 1 — 2n 



E n = (n + \)h - (n + h) (h 2 cu 2 /W m ) -V m lu = a(2V m /m)* 

where L b n , b (z) are generalized Laguerre polynomials. Experimental data show 
that these eigenvalues give very accurate values for the vibrational levels of 
nearly all molecules, but again the condition u n (0) = is not satisfied. 



If one considers that by substituting the potential (5) "we are left with the 
equation of a one-dimensional harmonic oscillator centered at r = r m " [3, Ay], no 
boundary condition is required at the origin, so that in particular perturbation 
theory for a one-dimensional harmonic oscillator is applied to the third-order 
term. 

Otherwise, arguments have been proposed to justify that the fact that u n (0) ^ 
is of no significance. According to Pauling et al. the solutions of the radial 
equation must be zero at r = and +oo, whereas the solutions of the harmonic 
oscillator must be zero at r = -oo and +oo, but because of the rapid decrease 
in the harmonic oscillator functions outside the classically permitted region, "it 
does not introduce a serious error to consider that the two sets of boundary 
conditions are practically equivalent" [4, See also 5,6,7]. As noted by Schiff, the 
eigenvalues of the radial equation are those of a linear harmonic oscillator "if 
the domain of r is extended to — oo" [8], which amounts to substituting a one- 
dimensional problem for a three-dimensional problem. Or else, one considers 
that, as (6) is the radial equation for zero angular momentum states, the exact 
solutions must be rigorously zero at r = 0, but the solutions (7) are "practically 
zero at the origin" [3,Fvn]- It should be noted that one often considers indeed 
that the exact solutions of (6) can be obtained only if the variable is taken from 
— oo to +oo, or if the solutions vanish at the origin. In fact, (7) are the exact 
solutions of (6), that is, the solutions in R 3 (r >0), as is easily checked. 

Similarly, the value of the Morse potential is sometimes considered as so 
large at r = that the Morse eigenfunctions are "effectively zero for r < 0" [9]. 
In fact they are non zero at r = 0, but according to Morse since "in every case 
r^ will be extremely small . . . this discrepancy will not affect the values of the 
energy levels" [2, See also 6,7]. 

But if the two boundary conditions, u n (0) = and u n (0) « 0, are practically 
equivalent, it merely means that the former condition is not required. First- 
order singular functions are not indeed "practically eigenfunctions" of H, they 
are eigenfunctions of H in R 3 /{0}, but not in R 3 . In other words, whether or 
not m(0) is very small, if it is non zero, ^/„ behaves at the origin like 1/r, and 
what these examples show is thus that energy levels are not affected because 
that they are not given by eigenfunctions of H. By changing the condition at the 
origin, we change indeed the operator. 

3 The energy levels and the operator Hj. 

We must thus determine the operator which has singular eigenfunctions of the 
form ^> n = \u n (r)/r\Y^ (6,ip) and corresponding eigenvalues E n , such that for 
£ = u n (r) and E n are solutions of (4) with u n (0) 7^ 0. Now, when ^ n and 
u n (r), which are of C°°(R 3 /{0}), are singular at the origin, they define in R 3 
distributions called pseudofunctions, denoted by Pf.\& n and Pf.u„(r) [10], and 
the latter are not a solution of (2) [1]. The radial equation which has in R 3 
singular solutions Pf.u ra (r) given by the root A = — I is the extension in R 3 of 



(2), it is the equation [1] 

h 2 d 2 £(£+i)h 2 



Pf. 



2m dr 2 



2mr 2 



+ V(r) 



u n (r) = E n Pi.u n (r) 



(10) 



Its solutions, given by both roots A = £ + 1 and A = —£, are either functions 
u n (r) which behave at the origin like r e+1 , or pseudofunctions Pf.«„ (r) which 
behave at the origin like Pf.r . (The symbol Pf. is used when we have either 
the former or the latter). When the solutions of (10) are given by the roots A 
— £ + 1 and A = (i.e. A = — £ for £ = 0), the symbol Pf. is useless and (10) is 
written in the form (2). Just as ^„ defines in R 3 the pseudofunction Pf.^n = 
Pi.[u n (r) /r]Y^ (9,(f), its Hamiltonien in R 3 /{0}, which is the Hamiltonian in 
the sense of the functions 



r 



tf_£_ £{£ + l)h 2 
2m dr 2 2mr 2 



V(r 



Un(r)Yf(e,<p) 



r> 



defines in R 3 the distribution 



Pf.tftf„ = Pf. i- 



2m dr 2 



£(£+l)h 2 
2mr 2 



V(r) 



u n (r)Yt (e,<p)\ (11) 



Let Hd the operator whose action on a pseudofunction Pf.'J'n amounts to 
the embedding in the distributions of R 3 of the Hamiltonian of ^ n in R 3 /{0}, 



tf d Pf.*„ = Pi.HV r 



(12) 



If Pi.u n (r) is a solution of (10), whether it is given by the root X = £ + I, 
or by the root A = —£. by substituting (10) in (12), we obtain the eigenvalue 
equation of Hd , 

HdPi.^a = E„Pf.*„ (13) 

where E n is the eigenvalue of (10) that belongs to Pf.M„ (r). 

The normalizable solutions are given by the roots A = £+ 1 W, and A = 0, 
that is, A € N. As they are less singular at the origin than 1/r 3 , the symbol Pf. 
can be dropped [10], and (13) is written 



flrftfr. 



XeN 



(14) 



In particular for £ — 0, ^ n — (4w) 1 / 2 [u„(r)/r], whether u n (0) = if A = 
1, or u n (0) ^ if A = 0, in both cases 



H, 



1 u n (r) 
^/47r r 



1 F u n (r) 



4tt 



A = 



(15) 



where u n (r) and E n are solutions of the radial equation (4). Hd is thus the 
required operator. 



In order to compare the operators Hd and H, let us recall that the Hamilto- 
nian of Pf.^„= Pf.r A_1 ^2 k>0 a k r k Yj? (9, ip) is the Hamiltonian in the sense of 
the distributions given by [T\ 



r 



h 2 d 2 e(£ + i)h 2 



2m dr 2 2mr 2 



V(r) 



2 



"„(>■)¥," (0. T -)\-^Q x ^{S) 



Qx,i (S) = ]T akX P (k)Be, m C p{k) r e Y$ (9, V )X>^5 p = - - - "* 



fe=0 



_/l ifpeN _ \-2t (4p+l)7r 3 / 2 

Xp ~\0 ifp^N ^ P ~4 P +1 Up ~ 22p-yr(p + 3/2) [W > 

where A p is the iterated Laplacian, <5 the Dirac mass, where the even, resp. 
odd, terms of the sum are zero , if A — £ is odd, resp. even, and where Q\t (5) 
^ if and only if there is at least one coefficient a k ^ for which k + A — I is 
an even negative integer [1]. As 

(H d - ff )Pf.tt„ = [Pf., H] tf „ = ^-Q XJ (5) (17) 

2m 

the operators H d and H differ when they act on pseudofunctions for which 

the operators Pf. and H do not commute, or for which Q\ t i (S)^ 0. Considering 

that energy levels are given in any case by the operator H amounts to identify 

the operators Hd and H, and hence to take no account of the noncommutation of 

the operators Pf. and H. As the operators Hd and H are equivalent in R /{0}, 

and in R when the Laplacian is the Laplacian in the sense of the functions, it 

comes from the confusion between the Laplacians in R /{0} and in R , or/and 

between the Laplacians in the sense of the functions and in the sense of the 

distributions [11]. By substituting (10) in (16), we obtain 

FPf.tf„ = £„Pf.*„ - ^Qxj (5) (18) 

2m 

When the solutions of the radial equation are given by the root A — £+ 1, 

p = - (k + l)/2 <^ N Vfc, so that Qi+\^ = 0, \I/„ behaves at the origin like 

r e Yf (9,(p), it is an eigenfunction of H. When the solutions of the radial are 

given by the root A = —£, as a ^ and p = £ e N for k = 0, Q-g.i ^ 0, Pf.^n 

behaves at the origin like Pf.r~^ +1 ) Y/ 4 {9, (p), it is not an eigenfunction of H [1]. 

As Qo.o = -vSu„(0)(5, in the case of normalizable solutions (18) is written 

H^ n = E n V n \ = £+l 

H lUr^ = l En Un(r)_ + ^Vl Un{0)5 A = Q (19) 



y/Air r ^/Ak 



Comparison between (14), (15) and (19) shows that energy levels, which 
are given by eigenf unctions of H when u n (0) = 0, are given in any case by 
eigen] 'unctions of Hd- 

4 The Hamiltonian, the Schrodinger operator and 
the operator H^. 

The differential form 

h 2 
H = - — A + V (20) 

2m 

is usually called Hamiltonian or Schrodinger operator. Now, a differential 
form such as (20) defines an operator only on a given class of functions, so 
that the same differential form can define different operators [See e.g. 12]. Any 
function \I/, or pseudofunction Pf.^, belongs to the domain T> on which the 
Hamiltonian, denoted by H, is defined by (20). On the other hand, the operator 
involved in the Schrodinger equation, or Schrodinger operator, denoted by H s , 
can be defined on the set V s of the functions which are any superpositions 
of eigenfunctions of H. As T> s is a subset of T>, and as the two operators are 
equivalent on T> s , that is, H^> = H s $> if \& e T> s , the Schrodinger operator is a 
restriction of the Hamiltonian. As long as one confines oneself to solutions of the 
Schrodinger equation, these two operators, which have the same eigenfunctions, 
are equivalent, so that the distinction is not necessary, and in fact it is not made. 
Nevertheless the distinction can be made with respect to the operator Hd- 

The latter is defined on the same domain T> as H. Moreover (13) and (18) 
show that the eigenfunctions of H, or of H s , form a subset of the eigenfunctions 
of Hd, so that if \I> G T> s , then Hd^ = H s ^ . The operator Hd is thus also an 
extension of the Schrodinger operator. But the fact that the eigenvalue equation 
of H d is an extension of the eigenvalue equation of H to the pseudofunctions 
such that <5a,£ {$) 7^ does not mean that H d is an extension of H. Eq.(17) 
shows indeed that they are different when Q\^i (8) ^ 0, or when the operators 
Pf. and H do not commute. It follows that the operators H and H d are two 
different extensions of the Schrodinger operator. 

The condition for the operator H (and then H s ) to be self- adjoint is usually 
written u(Q) = 0. However, according to Merzbacher the condition that the 
Hamiltonian must be self- adjoint implies that any two physically admissible 
eigensolutions of (2) must satisfy the condition [13] 

lim_o(^-^)=0 (21) 

This condition is satisfied, either for w(0) = and u '(0)^ 0, in which case 
$ is an eigenfunction of Hd and of H, or for u '(0) = and u(0) ^ 0, in which 
case ^ is an eigenfunction of Hd, but not of H It follows that (21) is in fact the 
condition for Hd to be self-adjoint. 



The operator Hd can also be defined, by using the linearity of the operator 
Pf., from the operator Pf.A introduced in [11]. As 

h 2 
Pf.H = Pf.( A + V) 

2m 

by similarly defining the operator A^ as 

A d Pf.* = Pf.A* (22) 

we have 

h 2 

H d = -—A d + V 

2m 

5 A physically not required condition. 

The boundary condition w„(0) = is definitely stated as a condition required 
to obtain energy levels because it is the condition for \I/„ to be an eigenfunction 
of the Hamiltonian, but in practice this condition is not respected. It can 
be imposed when the solutions are determined with the help of numerical or 
approximation methods, but not when energy levels are determined with the 
help of analytic solutions. 

In the case of diatomic molecules, in the parabolic like in the Morse potential, 
the fact that u n (0) ^ is regarded as of no significance owing to the fact that 
the origin being some distance outside of the classically allowed region, m„(0) is 
"very small". As the parabolic potential (5) is an expansion of the potential in 
the neighborhood of the minimum by neglecting the terms of order higher than 
2, it remains close to the physical potential only in the neighborhood of the 
minimum. The Morse potential (9) is closer to the physical potential, and in a 
region more extended on both sides of the minimum. Comparison between the 
measured vibrational levels £ n and the eigenvalues E n of the radial equations 
involving these potentials shows that the agreement between these eigenvalues 
and the energy levels decreases with n, and that the number of eigenvalues which 
are close to energy levels is all the greater as the region in which V(r) rj V(r) 
is more extended. It means that this agreement essentially depends on u n (r) in 
the region where V(r) is close to V(r), where one therefore expects that u n (r) is 
close to w n (r), the region where nearly the whole wave function is concentrated. 
It follows that if energy levels are not affected by the fact that m„(0) ^ 0, it 
does not come from the fact that it„(0) is "very small", but from the fact that 
the origin is outside of the region where V(r) w V(r), where u n (r) w w n (r), 
and where the wave function is concentrated, that is, in a region where u n (r) 
has little or none influence on the fact that E n is close or not to £ n . 

Besides, in the case of the potential (5) u n (r) always vanishes inside of the 
classically allowed region. It means in particular that the condition u n (0) = 
can be satisfied only inside of this region. 



Consider a physical potential V {r), and a supplementary potential, or a 
perturbation, W (r) , for which one substitutes the spherical oscillator V (r) = 
\muJ 2 r 2 and the linear potential W(r) = - Cr, so that the approximation of 
the physical potential V (r) = V (r) + W (r) is the potential V(r) = V (r) + 
Fr(r) of the form (5), that is, 

V (r) = Jmw r — Cr = \mui (r — r m ) + V^j 

r m = C/muj 2 V m = C 2 /2muj 2 (23) 

As C and then r m have any values, the origin is or not outside of the 
classically permitted region. When fir m coincides with a zero of a Hermite 
polynomial of order N. that is, Hn {(ir m ) = 0, then ujy (0) = 0, but u n (0) ^ 
for n < N. Insofar as the agreement between the eigenvalues of the radial 
equation and the energy levels decreases with n, the eigenvalues E n with n < N 
that belong to solutions u n (r) which do not satisfy the condition u n (0) = must 
be closer to an energy level than En that belongs to a solution ujy (f) which 
satisfies this condition. The latter then cannot be considered as the condition 
required to obtain accurate, or better estimates of energy levels. 

In the absence of supplementary potential, or of perturbation, the approxi- 
mation of the physical potential is the spherical oscillator, V(r) ~ V(r), with 

V (r) = \muj r 



u n (r) = N n exp(-^ 2 r 2 )H n (^r) E n = (n + \)hu 

u n (0) = N n H n (0) 

As Hermite polynomials satisfy the condition H2 P (0) ^ 0, i?2p+i (0) = 0, 
half the square-integrable solutions of the radial equation satisfy the condition 
u n (0) = 0, and hence half the energy levels are given by eigenvalues of H. If the 
supplementary potential, or the perturbation, is a ID potential W(x) for which 
one substitutes the linear potential W(x) = - Cx, the potential V(r) = V (r) 
+ W(x) used as an approximation of the physical potential V(r) = V (r) + 
W(x) is still a spherical oscillator, that is, 

V (p) = \mijj 2 r 2 — Cx — \rwjj 2 p 2 — V m 

p=[(x- C/mio 2 ) 2 + y 2 + z 2 ] V m = C 2 /2mio 2 

Half the energy levels are still given by eigenvalues of H. But when the 
supplementary potential, or the perturbation, is the 3D potential W(r) for 
which one substitutes the linear potential W(r) = - Cr, the potential used 
as an approximation of the physical potential is the potential (23), so that 
the solutions of the radial equation, which are given by (7), do not satisfy the 



condition «„(0) = 0, and then that the energy levels are not given by eigenvalues 
of H, as it is the case for diatomic molecules. It shows that whether or not energy 
levels are given by eigenvalues of H only depends on the potential used as an 
approximation of the physical potential - and hence on the latter. 

Perturbation methods, in particular variational methods, apply to the eigen- 
functions of H on which the theory is based. By substituting Hd for H these 
methods are extended to the eigenfunctions of Hd which behave at the origin 
like 1/r. It is in agreement with the fact that energy levels are obtained with 
the help of such functions. Similarly, situations which are described with the 
help of functions which behave at the origin like 1/r, require, to hold in R 3 , 
and not only in R 3 /{0}), that one substitutes Hd for H. This is the case of 
ingoing, outgoing and standing spherical waves which are described with the 
help of Green functions. As the latter behave at the origin like 1/r, they are 
eigenfunctions of H in R 3 /{0}), but not in R 3 , so that this description does 
not hold at the origin, which is physically unsatisfactory. As Green functions 
are eigenfunctions of Hd in R 3 , by substituting this operator, the description of 
spherical waves in terms of Green functions holds in the whole space. 

Moreover, wave functions divergent at the origin, but normalizable, are en- 
countered in the relativistic theory of the hydrogen atom. As the relativistic 
generalization of the Schrodinger equation is the Klein-Gordon equation which 
governs the wave function in the absence of electromagnetic field for a relativis- 
tic particle with spin zero, if the Laplacian must be taken in the sense of the 
distributions in this equation, one will be similarly led to substitute for A the 
operator Ad defined in (22). 

6 Conclusion. 

One usually considers, by putting forward different arguments, that the reso- 
lution of the Schrodinger equation must be supplemented with the boundary 
condition u n (0) = 0, but when this condition is not satisfied, one admits, by 
putting forward other arguments, that this is of no significance. In fact the reso- 
lution of the Schrodinger equation need not be supplemented with any boundary 
condition [1], but the fact that the condition w„(0) = is not satisfied is sig- 
nificant, since it means, whether or not u n (0) is very small, that energy levels 
are not given by solutions of the Schrodinger equation. As energy levels are ob- 
tained by substituting theoretical and then approximate potentials, it should be 
taken for granted that in approximate solutions the boundary condition u n (0) = 
0, and hence the operator H, or H s , are not required to obtain energy levels, and 
hence that the operator Hd, less stringent, is sufficient. The point is that the 
only potentials that are exactly known, and then used to obtain energy levels, 
are theoretical and then approximate potentials. 
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